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ABSTRACT 

An  IBM  7090  computer  program  is  described  which  simulates  vehicle  assign 
ments  to  a  priority  ordered  sequence  of  cargo  units.  The  schedule  thereby 
generated  is  used  to  assess  the  lift  capability  of  an  arbitrary  vehicle  inventory 
or  to  compare  one  vehicle  inventory  with  another.  The  model  can  thus  be  used 
to  determine  the  logistic  feasibility  of  an  operations  plan;  and  if  infeasibility  is 
demonstrated,  indicate  where  logistic  augmentation  or  modification  is  most  de¬ 
sirable.  Flowcharts,  a  listing  of  the  FORTRAN  program,  and  sample  inputs 
and  outputs  are  included. 
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I.  INTRODUCTION 


The  Logistic  Feasibility  Scheduling  Model  (LFSM),  as  presently  formulated 
can  be  utilized  with  either  of  two  functions  in  mind.  These  can  best  be  described 
as  requirements -oriented  or  capabilities -oriented .  In  the  former  application,  the 
leasibility  test  accepts  the  schedule  of  requirements.  It  then  assesses  the 
capacity  of  specified  logistic  resources  to  meet  the  schedule.  Requirements  that 
cannot  be  satisfied  are  noted  as  indicating  the  need  for  augmentation  of  capabilities. 
If  capabilities -oriented,  the  feasibility  test  considers  the  logistic  resources  inputs 
as  binding,  but  construes  the  requirements  schedule  to  be  indicative  of  a  desired 
time  sequencing  rather  than  as  a  firm  schedule. 

The  LFSM  is  designed  to  consider  all  cargo  units  as  comprising  a  single 
theater-wide  priority  list.  Total  transportation  resources,  both  organic  to  the 
theater  of  operations  and  extra -theater,  are  surveyed  in  the  determination  of  the 
lift  schedule. 


II.  GENERAL  DESCRIPTION 

The  LFSM simulates  the  following  situation:  A  given  set  of  cargo  units 
must  be  transported  from  their  origin  ports  to  their  delivery  ports  in  the  shortest 
time  by  a  given  set  of  vehicles.  Each  cargo  unit  has  associated  with  it  an  origin 
port,  a  delivery  port,  a  delivery  time  limit  (DTL),  a  tonnage,  a  priority,  and  a 
time  of  availability.  For  each  day  in  the  total  time  period,  every  cargo  is  con¬ 
sidered  in  turn  (cargoes  are  not  considered  until  the  time  they  become  available) 
on  a  priority  basis?  Once  the  highest  priority  cargo  available  is  chosen,  all 
compatible  vehicles  are  considered  to  determine  the  earliest  time  that  each  can 
make  the  lift  and  be  accepted  at  the  delivery  port.  If  the  earliest  of  these  times 
is  less  than  or  equal  to  the  required  delivery  time  limit,  a  movement  will  be 
scheduled;  otherwise,  not.  Normally  the  vehicle  with  earliest  delivery  time  will 
be  unique  and  is  automatically  scheduled.  If  more  than  one  vehicle  has  the  same 
earliest  delivery  time,  within  the  DTL,  the  one  with  greatest  cargo  capacity 
is  scheduled. 

If  one  vehicle  does  not  suffice  to  move  the  whole  cargo  unit,  the  weight  of 
the  unit  is  reduced  by  the  capacity  of  the  vehicle  assigned,  and  the  above  cycle  is 
repeated.  In  this  way  it  is  possible  for  a  fraction  of  a  cargo  unit  to  be  scheduled 
even  though  delivery  of  the  full  unit  is  not  feasible  within  the  required  delivery 
time  limit.  On  the  other  hand,  if  a  cargo  unit  or  its  remainder  does  not  fill  the 
assigned  vehicle  and  there  is  excess  capacity,  other  cargo  units  at  the  same 
origin  port  are  considered  in  order  of  priority.  If  one  of  these  cargo  units  has 
the  same  delivery  port  as  the  first  unit  assigned,  is  compatible  with  the  assigned 
vehicle,  is  available  by  the  time  the  vehicle  arrives  at  the  origin port,will  fit 
completely  within  the  remaining  vehicle  capac  ity,  and  can  be  delivered  within 
its  DTL,  it  is  assigned.  Such  a  unit  may  be  well  down  the  priority  list  as  com¬ 
pared  to  the  first  unit  assigned,  but  its  transport  does  not  delay  any  higher 
priority  movement.  Fractional  units  are  not  assigned  in  these  circumstances, 

♦In  this  model  both  the  day  of  availability  and  the  priority  of  the  cargoes  are  used  to 
determine  which  cargo  is  chosen.  There  exists  another  model,  computer  program 
38  -64P1,  in  which  just  the  priority  of  the  cargoes  are  used  todetermine  which  cargo 
is  chosen. 
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to  prevent  a  single  cargo  from  being  delivered  at  two  widely  separated  points  in 
time,  since  the  model  will  designate  for  the  next  lift  the  cargo  unit  of  highest 
remaining  priority  whether  or  not  it  is  the  partially  delivered  cargo  unit. 

The  LFSM  thus  generates  a  feasible  schedule  of  delivery  for  priority  ordered 
cargo  units,  transportation  being  scheduled  only  for  that  fraction  of  each  unit 
which  can  be  delivered  by  the  stipulated  delivery  time  limit.  Cargo  units  (partial 
or  complete)  that  cannot  meet  the  desired  delivery  schedules  are  regarded  as  im¬ 
posing  requirements  for  vehicle  augmentation  and/or  port  improvements. 

The  LFSM  may  be  used  to  assess  lift  feasibility  in  any  of  the  three  modes; 
airlift,  sealift,  or  combined  lift.  In  the  latter  case,  provision  is  made  by  the 
model  for  predetermined  assignment,  by  the  logistics  planner,  of  those  specific 
cargo  units  to  be  transported  strictly  by  air  and  those  strictly  by  sea.  This  is 
done,  by  use  of  a  cargo  compatibility  index.  (See  section  II  -  Inputs  to  the 
Program . )  Predetermined  cargo  unit  assignment  to  any  of  certain  specific 
vehicles  is  done  in  the  same  manner. 

The  outputs  from  the  computer  model  include  trip  index  quantities:  i.e., 
for  each  trip,  the  cargo  unit,  cargo  priority,  the  vehicle  identification  number, 
origin  port,  departure  time,  delivery  port  (destination),  arrival  time,  and  ton¬ 
nage  for  each  cargo  carried .  Other  outputs  include  the  tonnage  of  each  cargo 
unit  remaining  undelivered  along  with  its  DTL;  a  matrix  of  cumulative  unused 
vehicle  capacities  for  each  type  of  vehicle  at  each  delivery  port;  total  tonnage 
delivered  for  each  day  in  the  total  time  period;  and  final  delivery  times  for 
each  cargo . 

III.  INPUTS  TO  THE  PROGRAM 


Address 

1-700 

Symbol 

CP. 

l 

Description 

Tonnage  of  i^1  cargo  unit  (l^i  s700) 

701-1400 

PC. 

l 

Priority  of  i^  cargo  (PC  =  1  for  highest, 

1401-2100 

YP. 

l 

PC.  =  2  for  next  lower)* 

Day**  of  availability  of  iL  cargo 
(0  <  YP.  <;  NDAYS)  Note  YP.  t  0 

*  It  is  possible  for  two  or  more  cargoes  to  have  the  same  priority  number.  If 
this  is  the  case,  the  computer  program  will  consider  the  one  with  the  lowest  address 
as  the  highest  priority. 

**  All  times  are  entered  in  days  and  fractions  of  a  day  according  to  the  following 
convention:  0  <  DAY  s  1  is  considered  as  first  day 
1  <  DAY  s  2  is  considered  as  second  day 
NDAYS  -  1  <  DAY  s  NDAYS  is  considered  as  last  day 
i.e. ,  if  DAY  =3.7,  we  are  seven -tenths  of  the  way  through  the  fourth  day. 


III.  INPUTS  TO  THE  PROGRAM  (cont’d) 
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\  * 

Address 

Symbol 

Description 

• 

1 

2101-2800 

ZP. 

l 

Delivery  time  limit  (DTL)  for  delivery  of 
itB  cargo  (days) 

l  ■ 

2801-3500 

OP. 

i 

Origin  port  of  i^  cargo 

3501-4200 

DP. 

l 

t  h 

Delivery  port  of  l  cargo 

4201-4900 

TC. 

l 

Type  of  i1"*1  cargo  (1  £TC^15) 

1 

4901-4915 

NAMAX 

P 

Maximum  number  of  air  slots  per  day  at 
ptB  port  (1  <p  £  15) 

' 

.  4916-4930 

/ 

NSMAX 

P 

Maximum  number  of  ship  slots  per  day  at- 
th 

p  port 

i 

4931 

NCARGS 

Number  of  cargo  units  (1  £NCARGS  £700) 

4932 

NVEHS 

Number  of  vehicles  (1  £  NVEHS  £500) 

1 

4933 

NDAYS 

Number  of  days  in  the  total  time  period 
(1  £  NDAYS  £90) 

1 

* 

4934 

5001 -5500 

* 

ID 

Vk 

Identification  numberof  this  computer  run(ID 
is  an  integer  such  that  1  £  ID  £  32767) 

ID  number  of  k1*1  vehicle  (1  £k  £500) 

(for  aircraft,  Vk  £  1000,  for  ships 

5501-6000 

VT, 

Vk  > 1000) 

Type  of  kth  vehicle  ( 1  £  VTk  £  1 5) 

I 

6001-6500 

wvk 

th 

Day  of  availability  of  k  vehicle 
(0<WVk  £  NDAYS)  Note  WVfc  f  0 

i 

| 

1 

6501-7000 

BVk 

Port  of  availability  of  ktB  vehicle 
(1  £  BVk  £  1 5) 

fU 

1,  1,  33-15,15,33  IV  Compatibility  of  v  type  of  vehicle  with 

v’  c  th 

c  type  of  cargo  (if  compatible,  IV  =  1, 

v ,  c 


if  not,  IV  =0) 

V)  ^ 

1, 1, 34-15, 15, 34  IBV  Compatibility  of  vtB  type  of  vehicle  with 

V»P  th 

p  port  (if  compatible,  IBVv  n  =  1,  if  not 

mvy>p  =  °) 
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III.  INPUTS  TO  THE  PROGRAM  (cont'd) 


Address 

Symbol 

1 ,  1,35-15, 15, 35 

DELV 

v,P 

] ,  1 , 36-15, 15. 50 

EV  * 

vfp.q 

1 , 1,51-15,  15, 65 

HV 

v,p,q 

Description 

th 

Processing  time  for  v  type  of  vehicle 
th 

at  p  port  (days) 

Transit  time  for  v  type  of  vehicle  from 

pth  port  to  qfl1  port  (days) 

Capacity  of  v1*1  type  of  vehicle  from  p^ 
th 

port  to  q  port  (tons) 


IV.  USER'S  INSTRUCTIONS 


See  the  list  of  input  parameters  on  page  2.  The  user  should  submit  a  data 
submittal  form  (see  appendix  C).  These  forms  have  space  for  the  submitter's 
name,  the  date,  the  program  number  (38-64P),  the  security  classification  if 
any,  and  the  parameter  addresses  and  values. 


Input  flexibility  has  been  attained  by  allowing  the  user  to  vary  any  or  all  of 
the  parameters  in  a  computer  run.  In  any  one  set  of  data  the  parameters  re¬ 
main  fixed,  of  course,  but  there  is  no  programmed  limit  to  the  number  of  data 
sets  a  user  may  submit  in  a  run.  The  only  restriction  is  that  each  data  set 
must  terminate  with  one  blank  card,  and  the  last  set  in  the  run  must  terminate 
with  two  blank  cards, 

A  further  advantage  enjoyed  by  the  user  is  that  for  each  data  set  after  the 
first,  he  need  submit  only  those  parameter  values  in  a  set  that  are  different 
from  those  in  the  previous  set.  This  is  accomplished  by  identifying  each 
parameter  by  a  unique  "address"  in  the  computer  memory  (see  appendix  E). 
Initially  every  address  is  cleared  to  zero  so  that  only  non -zero  input  parameters 
need  be  entered. 


V.  SAMPLE  PROBLEM 

For  illustrative  purposes,  the  following  sample  situation  was  fabricated: 

One  large  cargo  unit  must  be  lifted  from  its  origin  port  to  a  delivery  port  by  a 
group  of  four  vehicles  -  two  ships  and  two  aircraft.  The  specific  parameters  of 


*  As  an  example  of  how  the  subscripts  are  related  to  the  data  addresses:  If 
the  transit  time  of  the  second  type  of  vehicle  from  the  fourth  port  to  the  third 
port  is  5  days,  then  EV0  .  ,  =  5  and  address  2,4,38  has  the  value  5.  (See 

section  IV  -  User’s  Instructions  .) 


the  problem  are  shown  in  the  problem  submittal  form  in  appendix  C.  The  out¬ 
put  is  presented  in  appendix  D. 

VI.  KEYPUNCH  INSTRUCTIONS 

See  the  list  of  input  parameters  on  page  2  and  the  sample  problem  submittal 
form  in  appendix  C.  Input  card  formats  are  described  in  appendix  E. 

VII.  OPERATOR'S  INSTRUCTIONS 

Run  under  control  of  the  Bell  System  on  the  IBM  7090.  No  special  instructions . 
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CLEAR  MATRIX  OF  UNUSED  VEHICLE  CAPACITIES  I  J  SAVE  BV  for  k  =  1,2.3 .  SOU  IN  A NOTH  I  R  ARRAY  CALLED  AV.  |  IMTIALI7E  TRIP 


LEGIST  1C  FfASIHILITY  SCHtDOUNb  ME  l)F  l  38-  64P 

CIMtNSIEN  44(14630),  HV  (  15 ,  I  5 , 1 5  )  ,  li  V  (I  5 ,  t  5  ,  I  S  )  ,  Ofc  L  V  (  I  5  ,  I  5 )  , 

1  r- 1  l-i  V  (  lb,  lb)  ,  FIV{15,I5),  BY(-jQO),  *9(500),  VT(50O),  V(500), 
2FNSMAXI 15),  FNAMAX(ls),  TC  C  700 )  ,  L)P(70G),  ,;P(700),  ZP1700),  YPI700 
3),  PCI  700  ) ,  CPI  700  ) ,  AKU19C),  IP!  700),  UC(I5,15),  CM  700),  PK(7C0) 
4,  FNAI  15, 9J),  FNS(I5,9G)  , A9 ( 500 ) , X 9 ( 500 )  .ABCI15) 
fcUGI VALtNCc  (  C  <v  <  I  ) .  CPIIM,  ( 44 ( 70  I  )  ,  PCIIII,  (00(1401),  YP(I)), 
KCU2101),  ZPim,  (00(3501),  DP  I  1)1  ,  (00(4201),  TCI  1  )  ) 

2,  (  44  (  490  I )  ,  FNAMAX  (  1  )  ) ,  (00(4916),  FNSMAX(U),  (4U50C1),  9(1)), 
3(04(550  I  )  ,  VTIll),  i 44 ( 6C0  I ) ,  WV  (  1 )  ) ,  ( OOl  650  1  ) ,  RV(l)),  (<Jl'(7?01), 
4FIVID),  (0417426),  FlHV(l)),  (  04(  7651),  DtLVID),  (00(  7876),  fcVII 
5)),  (44(  11251),  HVI  1 )), (40(2801 ) ,0P(  1  )  ) 

C£  5  I  C  T  =  I,  15 
ARC ( IDT ) * ( ♦  1H  ) 

ABU  3 )  =  ( ♦ I FP  ) 

ABU  6  )  =  (  ♦  I  F  B  ) 

ABU  9 )  =  ( ♦ I  HR  ) 

ABU  12  )  =  (  ♦  I  FiT  > 

OB  10  1  1  1  =  1,  14630 
44(111 )=0.U 

CALL  DATA  (44,  15, 15,65,  1  ND > 

IF  (INU)  40,40,990 
IBSWCH=4Q(4935 ) 

IREBBT  =  0 
NCARGS=4U (493  1  ) 

NV£FtS  =  40  (  49  32  1 
NCAYS=44( 4933) 

NRUN=44  (  4934 ) 

HANC=( +6H38-64P) 

HAKD=( ♦  IH1) 

IF  (IBShCH)  50,50,60 
JCB Y=0 
MAM  =  0 
NSB  =  0 

JCBY= JCEY+ I 
MASS  I V=  3 

PRIM  52,HANL,NKUN, JCBY 

FORMAT  (  1)11,  39  X ,  18HINPM5  FOR  PROGRAM,  IX,  A6,  1  OX  ,  4HCASE  ,  I  6 , 25X  , 
14HPAGE,  1 6/ / 1 9  X ,  1H1,  9X,  1H2,  9X,  1H3,  9X,  1H4,  9X,  1H5,  VX,  IH6, 
29X ,  1H7,  9X ,  lhe,9X,lH9,ex,2HlO) 

NUM=M AM 
NSk=MJM 
MGM=  M  AM  ■»  1 
M  AM=MBM+9 

IF  (14630-MAM)  70,54,54 

ns  o=nsb ♦ 1 

IF  ( 4 4 ( N  S  0 ) )  5  7,55,  57 
IF  (NSE-MAM)  54,53,54 

PRINT  5  8,  NUM,  (4UILA),  LA=MOM,  MAM) 

FORMAT ( S X , I  5 , 2 X  r T OF IG.2 ) 

MASS  I V  =  MASS I V+  1 

IF  ( MASS  I  V- 55  )  53,  51,51 

HARD= ( ♦  IH  ) 

PRINT  61,  HAND,  NRUN 

FORMAT  (  1 H 1  ,  39X, 19HJUTPLTS  FOR  PROGRAM,  IX,  A6 , 1  OX , 4HC ASE , 16,25X1 
PRINT  80  , HARD 

FORMAT  (  A 1  ,  4  A  ,  4HTRIP,  5X ,  5HC  ARG.,  ,  5X ,  L’HPR  UR  I  T  Y ,  5X  ,  7H  VEH I  CLt ,  5X  , 
16H0KIGIN.5X  ,  UHDtPARTURt  TIMC  , 5X , 1  1HDL ST  I NA T I BN , 5X ,  1 2HARR I 9AL  TIM 
2C,5X, 7HT0NNAGE  //) 


i  ay 

90  PR  (  I  A  V  1  =PC (  JAY) 

CT (  I  AY  )  =CP (  1  AY ) 
fPl  I  AY )=0.O 
I  A Y= I  AY ♦  I 

IF  (  I  AY-NCAKGS)  VO, VO,  IOC 
100  1 0 A Y=  1 

NO  AKU  (  1  0  A  Y  )  =0 . 0 

I PORT  =  I 

120  FNAUPiiRT,  ICAY  )=FUAHAX(  IPORD 
FNS (  I  P0RT ,ILAY)=FNSMAX(  I P0K  T  ) 

I  P0R  T  =  1  PZRT  ♦  1 
IF  (  I  PISkT-  I  b  )  120,  120,  I  3C 
130  I0AY=ICAY*l 

IF  (  IOAY-VO )  110,110,140 
140  00  150  1= I, lb 

C0  lbO  J= 1, lb 
IbO  UC(I,J)=0.0 

02  160  K*1,500 
AV ( K  )  =HV ( K ) 

160  X  V  (  K  )  -  W  V  (  K ) 

JA  Y  =  0 
JC A Y=  1 

170  PR  I 2RT  =1.0 
100  PT  =  PRim 
I  A  Y=  1 

190  IF  ( YP(  IAYI-FL0A1FI  JUAY)  >  200,200,240 

200  IF  I PR( IAY1-PR10RT )  240,270,210 

210  IF  (Pr-PRI0KU  230,230,220 

220  IF  (PRIIAY1-PI)  230,240,241 

230  PT  =  PR (  I  AY ) 

240  IAY-IAY+1 

IF  (  I  AY-NCAKOS )  190,  I  >0,25C 
250  IF  (PT-PRIJRl )  070,670, 260 
260  PR  I 0R  T  =  PT 
G0  TO  180 
270  I Q=0P (  I  AY  I 
I  R  =  I)P  (  I  AY  ) 

Ifc  =  TCI  1  AY ) 

200  FVK  =  1000000 .0 

KAY-  1 

290  IC=VriKAY) 

IF  (  F  I  V  (  IC,  1 1 )  -  1 . 0  )  4  70,300,4  70 
30U  IF  CFIHVl IC, IQ1-1.0)  470,310,470 

310  IF  (FI0VI IC, IRI-l.O)  470,320,470 

320  I A=  AV ( K AY ) 

AVC=XV ( KAY ) +tV ( IC, IA, IC ) 

FVQ=MAX1F(YP( lAYI.AVO) 

330  IT=FVQ+.99V99 

IF  (  ir-NOAYS)  340 , 34  0 , 4  7  C 
340  NAY=FNA( 10, I T ) 

IF  (VIKAYI- 1000.0)  360,360,350 
350  NAY=FNS( 10, IT ) 

360  IF  (NAY)  370,370,300 

370  FVC=F VO ♦ 1 .0 

G0  T  0  330 

380  AVK=FVO  +  DFL  V ( lC,lU)+tV(IC,IQ,IR) 

390  IW  =  AVR+. 99999 


IF  (Ik-NCAYS)  400,400,470 
400  NAY=FNA( IR, I n ) 

II  ( V ( KAY )- ICOO.O)  420,420,4)0 
4  10  NA Y  =  F  NS ( IR,  I W) 

420  IF(NAY)  450,430,440 
430  AVR= AVR+ 1 .0 
Ok)  T 0  390 

-.0  IF  (AVR-FVK)  460 , 4  50 , 4  70 

4  jO  IF  (HVl IC,1C, IR)-HV(  IH, IC, IR) )  4/0,470,460 

460  I  M  =  K  AY 

I H  =  I C 
I  L  CJ—  1  T 
1  L  K=  I  W 

FI  VO-=FVC  +  DLL  V  (  IC,  Id 
FVR= AVR 
I  LA-  I  A 

470  KAY  =  K AY+  1 

IF  (KAY-NVCHS)  290,270,480 
480  IF  { FVR-ZPI  I  AY  )  )  490,490,240 

490  JAY=JAY+1 
IF AY=FVR 

IF  (FVR-FL0ATFI  IFAY)  !  495,495,491 

491  IF  (  FVR-I  FLkJATFl  IFAY)+0.CU)  492,495,495 

492  FFF=FL0ATF( IFAY)*0.0I 
00  T 0  496 

495  FFF=F VR 

496  AK  J  =  M  I N  I  F  ( C  T  ( IAY),  HV ( IH , IC , I R ) ) 

I  VEH  =  V  I IM) 

IPPPP  =  PKUkr 

IF  I  IR0B0r- I )  499,497,497 

497  PRINT  498, J AY, IVEH,  F  L  VQ  ,FI-f  ,  AK  J 

498  F0RMATI  4X , I  4 , 3 1 X  ,  I  4 , 23X ,  16 . 2,  2 7X  ,  F6 . 2 , 4X , F8. 2 ) 

G0  T«5  50  1 

499  PRINT  500,  J AY , I  AY  ,  I  PPPP , I VEH, I Q ,FL VO , I R , FFF , AK J 

500  F0RMATI  4X, 14.5X,  I5,dX, I a,8X, 14, 7X, 14, 12X.F6.2, 12X, 15, I0X,F6.2,4 
1X.F8.2) 

501  AKUI I LR  )  =  AKU ( ILRJ+AKJ 
TPI I  AY ) =F VR 

CT( IAY  )=CT<  IAY l-AKJ 
CAP=HV ( IH, IC, IR) 

C AP=CAP— AKJ 
IF  (CAP)  510,690,510 
510  SAY=PRI0RT 
520  ST=S AY 
MC=1 

530  IF  IPR(MC)-SAY)  570,600,540 

540  IF  (ST-SAY)  560,560,550 

550  IFIPR(MC)-ST)  560,570,570 
560  ST=PR(MC) 

570  MC-MC*  1 

IF  ( MC-NCAR GS  I  530,530,  580 
580  IF  (ST-SAY)  690,690,590 

590  SAY=ST 

00  T 0  520 

600  IF  ( CT ( MC ) )  570,570,610 

610  IF  (0P( VC)-FL0ATF(  10) )  570,620,570 

620  IF  (DP(8C)-FL0ATF(  IK) I  570,630,  570 

630  LE=TC(NC) 


B 


IF  (FIV(  IH,LE)-1.0)  570,  640,570 
IF  (CT(MC)-CAP)  650,650,570 
IF  (YP(MC  )  -FVC)  660,660,570 
IF  (FVR-ZP(MC)I  670,670,570 
AK J  =  CT ( MC  ) 

I  SSSS*SAY 

PRINT  680, MC, ISSSS.AKJ 
FORMAT  ( 13X,  I5,8X,I5,78X,Fd.2) 

AKU ( ILR ) =AKU( ILRUAKJ 
TP(MC) =FFF 
CT ( MC ) =0. 0 
PR ( MC ) -0 . 0 
CAP=C AP- AKJ 
IF  (CAP)  570,690,570 
IB=XV( IM)*. 99999 
LEAV=FLVQ*. 99999 
WDEL=FVR+OELV(  IH,  IR) 

IF  (VI IMI-1000.0)  700,700,710 

WDELSWDEL*2.0*EV(1H,ILA,IG)+2.0*EV(IH,IQ,IR) 

IW-WDEL+. 99999 

IF  ( I  LA- ID)  720,750,  720 

IF  (VI IMI-1000.0)  730,730,740 

FNAl ILA, I B ) =FN A( ILA, 18)- 1.0 

(.0  T 0  750 

FNS( ILA, IB) =  FNS  (  ILA, IB)- 1.0 
IF  (VC  IP)— 1000.0)  760,76C,770 
FNA( IG, ILU)=FNA( IQ, ILQ)-1.0 
G0  T 0  780 

FNS  (  10,  ILQ) =FNS( IG, ILQ)- 1.0 
IF  (LfcAV-ILQ)  800,800,790 
ILU= ILQ+  1 
G0  T 0  750 

IF  ( V ( I M ) -  1 000 . 0 )  810,810,820 
FNA( IR, ILR)=FNA(  IR, ILR)- 1.0 
G0  T 0  830 

FNSl IR, ILR)  =  FNS (  IR, ILR)- 1.0 
IF  (IW-ILR)  850,850,840 
ILR= ILR+  1 

IF  (ILK-90)  800,800,850 
AV( IM)=  IR 
XVI IM)=wDEL 

UC(  IH,  IR ) =UC ( IH, IRl+CAP 
IF  (CT (  I  AY )  )  860,860,851 
IRDBOT  =  1 
G0  TO  280 
P  R (  1  AY ) =0. 0 
IROBOT  =  0 
G0  T 0  240 
JOAY= JDAY+ 1 

IF  ( JDAY-NUAYS )  170,170,  880 
PRINT  890 , ( L , AKU I L 1 , L= 1 , NUAYS ) 

FORMAT  (  1 H I /45X , 29HT0NNAGE  DELIVERED  0N  EACH  CAY  III  SOX , 3HDAY ,  1  OX 
l,7hT0NNAGE//(50X, I3,9X,Fe.2)l 
PRINT  900 . I LC  .  TP ( LC ) , LC= 1 , NCARGS ) 

FORMAT  ( 1 H 1 /42X.35HFINAL  DELIVERY  TIMES  FOR  EACH  CARG0///5OX.5HCAR 
1  GO , 10X,4HTIME//(50X,I5,7X,F7.2) ) 

PR  I  NT  9  10 

F0RMAT (  IH 1/ 25X , 56HCARGC  TONNAGES  NOT  DELIVERED  WITH  REQUIRED  DEL1V 


6-6 


1ERY  TIME  ///3OX.5HCARG0, IOX, 7HT0NNAGE , 10X, I0HT1ME  LIMIT//) 

11  =  1 

IF  ( CT ( II))  960,960,940 
PRINT  950,1  I.CTI  I  I  ),ZP! II) 

FORMAT ( 30X , I5,9X,F9.2, I2X.F7.2) 

II*II  ♦) 

IF  I  I I-NCARGS)  930,930,970 
PRINT  980 

FORMAT  (  I M 1 /3SX, 54H  UNUSED  VEHICLE  CAPACITY  BY  VEHICLE  TYPE  FOR  EA 
1CH  PORT///  58X, 7HVEHICLE  // 1 8X,  1H),6X,IH2,6X, IH3.6X, I  HU, 

26X, 1H5.6X,  IH6.6X,  1H7.6X,  1H6, 6X, 1H9,  5X , 2H 1 0 , 5X , 2H1 1 , 5X , 2H1 2 , 5X , 2H I  3 
3,5X,2H14,5X ,2H15) 

00  989  L8= 1 ,  1b 

PRINT  981,  A8C(LB),LB,(UC(LA,LB),LA=1,15) 

FORMAT! IX, A 1,5X, I  4 , 4X , 14! F6 . 2 , 1 X ) , F  6 . 2 ) 

CONTINUE 
GO  TO  30 

PRINT  1000  , NRUN 

FORMAT  (40X , 18H***«*«END  OF  CASE  , 16, 6H**«**»/ 1H1 ) 

CALL  ENCJOB 
ENC 
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APPENDIX  D 


SAMPLE  PROBLEM  OUTPUT 


VtHICLE  ORIGIN  DEPARTURE  TIME  Ut  ST  I  NAT  I 2N  ARRIVAL  TIME  T  VNNAGt 


'iNAl'.L  DLL  t  VLfttU 

£,'*  u  AC  H  DAY 

DAY 

TONUAGfc 

1 

0. 

2 

0. 

1 

4.00 

4 

8.00 

5 

0. 

6 

60. OC 

( 

44.00 

8 

6. 00 

V 

0. 

10 

0. 

1 1 

4.00 

1? 

60.00 

11 

0. 

14 

48.00 

lb 

4.00 

16 

0. 

17 

0. 

IH 

60.00 

IV 

4.00 

20 

0. 

21 

8.00 

22 

40.00 

2.1 

4.00 

24 

60.00 

2S 

0. 

26 

0 . 

27 

12.00 

28 

0. 

2  V 

4  0.00 

10 

60.00 

11 

4.00 

52 

0. 

14 

0. 

14 

12.00 

lb 

0. 

16 

1  uO . 00 

17 

0. 

48 

4.00 

IV 

0. 

40 

ri.OO 

4  1 

0. 

4  2 

0. 

4  1 

0. 

44 

0. 

4  b 

0. 

46 

0. 

4  7 

0. 

48 

0. 

uv 

u. 

so 

0. 

si 

0 . 

s? 

0. 

SI 

0. 

S4 

0. 

Sb 

0. 

56 

0. 

5  7 

c. 

58 

0. 

SV 

c. 

60 

0. 

FINAL 

DELIVERY  TIMES 

FOR  EACH  CARG0 

CARGO 

TIME 

- 

1 

39.90 

CARG45  T0NNAGES  NUT 

DELIVERED  WITH 

REQUIRED  DELIVERY  TIME 

CARG0 

T0NNAGE 

TIME  LIMIT 

1 

344.00 

o 

o 

. 

o 

■S' 

tig. 


UNUSED  VEHICLE  CAPACITY  BY  VEHICLE  TYPE  F0R  EACH  P0RT 


o 

I 

UJ 

> 


*-000000000000000 


«••••••••  t  •  •••  1  • 

—OOOOOOOOOOOOOOO 


"n  1  •  1  *  •  •  •  *  »  •  •  t  1  1  • 

—000000000000000 


fSl  ••••••••«•••••  • 

—000000000003000 


—000000000000000 


o  ••••••••••••••  • 

—000000000000000 


000000000000000 


s  . . .  . . 

000000000000000 


000000000000000 


UJ 

t/* 

< 

#*•••♦•  t  #  I . . 

000000000000000 

LL 

<a 

o 

2 


000000000000000 


000000000000000 


000000000000000 


000000000000000 


000000000000000 


a 


*  ►- 
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APPENDIX  E 


DATA  SUBROUTINE 


1 .  Introduction: 

Many  computer  programs  require  the  flexibility  of  varying  any  or  all  of  the 
parameters  in  a  computer  run.  Although  FORTRAN  is  fairly  flexible  in  its 
arithmetic  and  control  statements,  its  input-output  statements  are  -quite  rigid. 

In  order  to  read  cards  for  instance,  considerable  effort  must  be  expended  by  the 
FORTRAN  programmer  in  writing  his  input  statements.  This  subroutine  eliminates 
some  of  that  tedium.  The  concept  of  a  "data  set"  is  used.  A  data  set  consists  of  a 
sequence  of  punched  cards  terminated  by  one  blank  card.  A  parameter  deck  for  a 
computer  run  may  consist  of  several  data  sets.  Such  a  parameter  deck  is  terminated 
by  two  blank  cards. 

2.  Parameter  Addresses: 

The  primary  advantage  of  this  subroutine  over  FORTRAN  input  statements  re¬ 
sults  from  the  use  of  "parameter  addresses."  An  address  is  a  relative  location  in 
the  computer  memory.  It  is  the  subscript  of  an  array  or  matrix.  For  example,  in 
an  array  called  X  ,  the  parameter  value  would  be  located  at  address  53  .  By 

using  the  parameter  addresses,  a  user  of  the  program  need  submit  only  those  param¬ 
eter  values  in  a  data  set  that  are  different  from  tr.ose  in  the  previous  set. 

Three  types  of  addresses  are  permitted  by  this  subroutine. 

(1)  A  numevic  address  consisting  of  one  to  five  characters,  each  of  which 
is  a  digit  0-9.  Such  an  address  (n)  refers  to  the  nth  element  in  a 
specified  array. 

(21  An  alpha  address  consisting  of  one  to  six  characters,  the  first  of  which 
must  be  alphabetic  (A-Z).  The  remaining  may  be  alphabetic  or  numeric 
(A-Z  or  0-9).  Such  an  address  refers  to  the  m*1  element  in  a  specified 
array  (1  i  n  i  26),  where  the  first  character  of  the  address  corresponds 
to  n  as  the  26  letters  of  the  alphabet  correspond  to  the  integers  1-26. 

(3)  A  matrix  address  consisting  of  two  or  more  numeric  fields  separated  by 
commas.  For  example,  the  address  53,  47  refers  to  the  element  ir  the 
53r<3  row  and  the  47“*  column  of  a  two-dimensional  matrix.  There  is  no 
limit  to  the  number  of  dimensions  in  a  matrix  address. 

3.  Input  Card  Format: 

A  standard  submittal  form  (see  attachment)  has  been  designed  for  the  analyst. 
This  form  provides  for  entering  parameter  values  with  their  associated  addresses. 
The  user  indicates  blank  cards  to  separate  data  sets.  The  keypunch  operator  has 
the  option  of  punching  one  address  and  value  per  card,  or,  if  the  addresses  are 
sequential,  of  punching  one  address  and  several  values  on  a  card. 


Only  columns  1-72  of  a  card  are  used.  Each  column  must  contain  one  of  the 
following,  a  digit  (0.-9),  a  "  +  ”  or  "  -  ”  sign  or  a  dash,  a  letter  (A-Z),  a  period, 
a  comma,  or  a  blank.  Each  punched  card  must  contain  one  parameter  address. 

The  address  may  start  in  column  1,  or,  if  desired,  may  start  in  a  later  column, 
provided  all  columns  before  it  are  blank.  The  address  is  terminated  by  at  least 
one  blank  column.  Only  one  address  is  permitted  on  the  card.  Succeeding  col¬ 
umns  contain  one  or  more  parameter  values,  each  separated  by  one  or  more 
blank  columns.  A  value  may  be  signed  or  unsigned.  The  length  of  the  value  field 
is  variable.  No  blanks  are  permitted  within  a  value  field.  A  value  may  be  punched 
with  or  without  an  exponent.  An  exponent  is  recognized  by  '  0  presence  of  a  plus 
or  minus  sign  (or  dash)  between  the  fractional  part  and  exponent  part  of  the  value. 
Decimal  points  (periods)  may  be  punched  in  either  the  fractional  or  exponent  parts 
of  a  value.  If  more  than  one  value  is  punched  on  a  card,  those  after  the  f  xst  will 
be  entered  at  sequential  addresses  relative  to  the  address  of  the  first  value. 

4.  Usage: 

A  data  set  is  read  by  the  use  of  the  statement: 

CALL  DATA  (X,  I) 

in  a  FORTRAN  program  for  the  IBM  7090.  The  argument  X  is  the  name  of  an 
array  in  the  program.  The  argument  I  is  an  indicator  set  by  the  subroutine. 

This  indicator  may  be  tested  by  the  main  program  upon  return  from  the  sub¬ 
routine.  It  will  have  a  value  of  0  or  1  or  2  . 

0:  The  subroutine  has  read  a  data  set.  The  main  program  will  normally 
proceed  to  operate  on  this  data . 

1:  The  subroutine  has  read  the  second  blank  card  which  terminates  the 
parameter  deck.  The  main  program  will  normally  terminate  at  this 
point. 

2:  The  subroutine  has  read  a  "bad"  data  card.  The  main  program  may 
terminate  the  run,  or  ignore  the  card  and  return  to  the  subroutine  to 
read  the  rest  of  the  data  set. 

If  the  cards  to  be  read  contain  matrix  addresses,  additional  arguments  must  be 
included  in  the  FORTRAN  calling  statement: 

CALL  DATA  (X,  Dr  D2>  D3 . D  .  I) 

where  Di  is  the  i1*1  dimension  of  the  matrix  X. 

5.  Method: 

See  the  attached  flow  chart.  DATA  reads  parameter  values  and  loading  addres¬ 
ses  from  cards.  If  sense  switch  5  is  up,  it  will  read  the  values  and  addresses  from 
tape  (unit  A2).  It  converts  the  values  to  floating  point  numbers,  and  stores  them  as 
elements  of  an  array  specified  in  the  calling  statement.  The  elements  are  specified 
by  the  addresses.  If  a  card  (or  tape  record)  is  read  which  contains  non-permitted 
characters  (see  input  card  format  above),  DATA  prints  tht  statement  "bad  data  card, " 
followed  by  an  image  of  the  card  itself. 
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SYMBOLIC  LISTING 


DATA 


A1 

A 


A2 


FAP 

entry  data 

SXA  XI  .  1 
SXA  X2.2 
SXA  X4.4 
1.4 


CAL 


ADD  CORE 
STO  XLOC 


AXT  1.1 
SXA  *+1.1 


CAL  **.4 

ana  mask 


tnz 

TXI 


STA 

AXT 

SXO 


STZ 

AXT 


*+2 

*-4.1.1 
SXA  EX1T.1 
TXI  *+1.1. -1 
SXA  *+1.1 
CLA  **.4 
STA  A 1 
STA  F  1 A 
112 
1.1 
**»1 
TSX  HHREAD.4 
PZE  CARD 
TRA  EXIT 
TRA  BAD 
STZ  ADDRES 
STZ  VALUE 
STZ  EXP 
POINT 
1.1 

SXA  FIELD.l 
AXT  13.1 
TNX  HH.1.1 
AXT  42.2 
SXA  COLUMN. 2 
LXA  COLUMN. 2 
TNX  A2.2.6 
SXA  COLUMN. 2 
LDQ  CARD+12.1 
RQL  36.2 
PXD  0.0 
LOL  6 

STO  CHARAC 
ORA  FLOAT 
FAD  FLOAT 
STO  NUMB 
AXT  42.4 
CLA  CHARAC 
CAS  TABLE+42.4 
TRA  *+2 
TRA  *+3 
T IX  *-3.4.1 


RETURN  INDICATOR  ■  1 
READ  A  CARD 


ADDRESS  *  0 
VALUE  «  0 
EXP  -  0 
POINT  »  0 
FIELD  ■  1 


COLUMN  6T  72 


COLUMN  «  COLUMN+1 
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mn 


Mr:  -.-'--o''' 


01 

?>'■  . ii  ■;  \  * 


. , 


TRA 

BAD 

LXA 

FIELD. 2 

TRA 

Fl+1.2 

TRA 

F7 

FIELO-7  (EXPONENT  FIELD) 

TRA 

F6 

FIELD-6  (VALUE  FIELD) 

TRA 

F  5 

FIELD-5  (BLANKS  AFTER  ADDRESS) 

TRA 

F4 

FIELD-4  (ARRAY  ADDRESS) 

TRA 

F  3 

FIELD-3  (ALPHA  ADDRESS) 

TRA 

F  2 

F 1  ELD-2  (NUMERIC  ADDRESS) 

FI 

TXH 

B  .  4  .4 1 

FIELD-1  (BLANKS  BEFORE  ADDRESS) 

F1A 

STZ 

*# 

RETURN  INDICATOR  ■  0 

AXT 

2.2 

FIELD  »  2 

SXA 

FIELD. 2 

TXH 

H.4.31  ' 

NUMERIC  CHARACTER 

TXH 

I  .4.28 

SIGN  OR  DASH 

TXL 

BAD. 4. 2 

AXT 

3.2 

ALPHA  CHARACTER.  FIELD  ■  3 

SXA 

f Feld. 2 

TXI 

•+1.4.-2 

ADDRESS  ■  NTH  ALPHA 

SXA 

ADORES. 4 

TRA 

B 

F2 

TXH 

J.4.41 

BLANK  CHARACTER 

TXH 

H.4.31 

NUMERIC  CHARACTER 

TXH 

K.4,28 

SIGN  OR  DASH 

TXH 

BAD. 4. 2 

TXH 

L.4.1 

PERIOD 

AXT 

4.2 

COMMA. FIELD  ■  4 

SXA 

FI.ELD.2 

AXT 

2.2 

SXA 

DIHENS. 2 

DIMENSION  ■  2 

CLA 

ADDRES 

TEST  ADDRESS 

TZE 

BAD 

TMI 

BAD 

F2A 

STZ 

ADDER 

ADDER-0 

TRA 

B 

F3 

TXH 

J.4.41 

BLANK  CHARACTER 

TXH 

B.4.31 

NUMERIC  CHARACTER 

TXH 

£.4.28 

SIGN  OR  DASH 

TXH 

B.4.2 

ALPHA  CHARACTER 

TXH 

L.4.1 

PERIOD 

TRA 

BAD 

F4 

TXH 

M  .4 .41 

BLANK  CHARACTER 

TXH 

N.4.31 

NUMERIC  CHARACTER 

TXH 

P.4,28 

SIGN  OR  DASH 

TXH 

BAD, 4, 2 

TXH 

Q'4,1 

PERIOD 

TRA 

T 

COMMA 

F5 

TXH 

a,  4. 41 

BLANK  CHARACTER 

AXT 

6,2 

FIELD  ■  6 

SXA 

FIELD. 2 

TXH 

U.4.31 

NUMERIC  CHARACTER 

TXH 

W ,4 ,28 

SIGN  OR  DASH 

TXH 

BAD, 4. 2 

TXH 

G.4,1 

PERIOD 

'•  -sTyyj 

'Vi &&  .v,.*?  • 
'  & 


..A--  ■  ..  ..-.  _  .  -  , 


'*•  ^  M 


i"  i  Jijup.  « ii.il 


I 


»1 


tra  bad 

TXH  X.4,41 
TXH  U  .  4  •  3  1 
TXH  Y.4,28 
TXH  BAD. 4. 2 
TXH  G  .4 . 1 
TRA  BAD 
TXH  2.4.41 
TXH  BB.4.31 
TXH  EE. 4. 28 
TXH  BAD .4.2 
TXL  BAD. 4.1 
AX  T  1.2 
SX A  POINT. 2 
TRA  B 

LDO  ADDRES 
MPT  H 1 0 
XCA 

ACL  CHARAC 
STO  ADDRES 
TRA  B 

TXH  B .4 .30 
CLA  ADDRES 
SSM 

STO  ADDRES 
TRA  B 
AXT  5.2 
SXA  FIELD. 2 
TRA  B 

TXH  LI. 4. 30 
CLA  VALUE 
SSM 

STO  VALUE 
TRA  LI 
AXT  1.2 
SXA  POINT. 2 
AXT  6.2 
SXA  FIELD. 2 
TRA  B 
AXT  5.2 
SXA  FIELD. 2 
TRA  S 
LDO  ADDER 
WPY  H10 
STQ  ADDER 
TSX  T 1 . 4 
MPT  CHARAC 
XCA 

ADD  ADDER 
STO  ADDER 
TRA  B 
TXH  R.4.30 
CLA  VALUE 
SSM 


BLANK  CHARACTER 
NUMERIC  CHARACTER 
SIGN  OR  DASH 

PERIOD 

BLANK.  CHARACTER 
NUMERIC  CHARACTER 
SIGN  OR  D*  SH 


PERIOD.  POINT  ■  1 


ADDRESS  ■  10  X  ADDRESS  +  N 


+  SIGN 

SET  SIGN  OF  ADDRESS 


FIELD  -  5 


♦  SIGN 

SET  SIGN  OF  VALUE 


POINT  •  1 
FIELD  ■  6 

FIELD  ■  3 

ADDER  ■  10  X  ADDER  ♦  N  X  PROD 


POINT  •  1 


STO  VALUE 
TRA  R 
AXT  1.2 
SXA  POINT. 2 
AXT  6*2 
SXA  FIELD, 2 
LXA  EXIT, 2 
TX 1  •♦1.2. “3 
PXA  0.2 
SUB  DIHENS 
TNZ  BAD 
TSX  T 1 .4 
CLA  ADDER 
SUB  PROD 
STO  ADDER 
TZE  BAD 
TMI  BM) 

ADD  ADDRES 
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